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Abstract. We consider the three-dimensional semi-relativistic Hartree model 
for fast quantum mechanical particles moving in a self-consistent field. Un- 
der appropriate assumptions on the initial density matrix as a (fully) mixed 
quantum state we prove, using Wigner transformation techniques, that its clas- 
sical limit yields the well known relativistic Vlasov-Poisson system. The result 
holds for the case of attractive and repulsive mean-field interaction, with an 
additional size constraint in the attractive case. 



1. Introduction and main result 

In this paper we aim to establish the classical limit as e — -> 0+ of the semi- 
relativistic Hartree system (or Schrodinger-Poisson system), i.e. 



(1.1) 



ied t p e {t) =[H e , P E (t)}, xeR 3 ,te 
-kAV s =rf(t,x), k = ±1, 



subject to an initial data p £ (0) = pg, where the reason for choosing k = ±1 respec- 
tively will be explained below and the Hamiltonian operator H £ is given by 



H £ := \/-e 2 A + l + V s (t, x). 

Here the pseudo-differential operator for the kinetic energy is simply defined via 
multiplication in Fourier space with the symbol ^/|e^| 2 + 1, for the semi-classical 
parameter < e *C 1, a dimensionless scaled Planck's constant (all other physical 
constants are rescaled to be equal to 1). This scalar pseudo-differential opera- 
tor is frequently used in relativistic quantum mechanical models as a convenient 
replacement of the full (matrix-valued) Dirac operator. In (|1.1|) we denote by 
p £ (t) £ Si(L 2 (R 3 )) the density matrix operator of the system, i.e. a positive self- 
adjoint trace class operator acting on L 2 (R 3 ). The particle density n £ (t, •) £ L 1 (M 3 ) 
is then obtained by evaluating the corresponding kernel p £ (t, x, y), which, by abuse 
of notation is denoted by the same symbol as the density operator, on its diagonal 
[Hj, i.e. n £ (t,x) = p £ (t 7 x,x). The system (jl.ip describes the mean field dynam- 
ics of (relativistic) quantum mechanical particles in a mixed state. Denoting by 
{V'fljeN a complete orthonormal basis of L 2 (R 3 ) we can decompose the kernel of 
p £ (t) £ 6i(L 2 (R 3 )) via 

p £ (t,x,y) = VM(t,i)I(t,j), 
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with (constant in time) Xj E i 1 , Xj > [TH]. Using this representation, we arrive 
at an equivalent system of countable many nonlinear Schrodinger-type equations 



(1.2) 



j ledt^ = y/-e 2 A + l ^ + V s (t, x)i% 
\ -nAV e = n £ (t,x), 
where the density is now given by 

The system (|1.2|) can be interpreted as the time-dependent model associated to the 
semi-relativistic Hartree energy 



£H(*)=£kin(*)+£pot(*) 



pot \ 



3& 

where for the second line we have used the three dimensional Green's function 
representation of the potential, i.e. 

V e (t,x) = -^—*n s {t,x). 
47r|a;| 

For future references we recall that the kinetic and the self-consistent potential- 
energy can also be shortly written in terms of density matrices as 

=tr(V- e 2A + lp e (i)), £ £ pot (t) = f tr(V E (t,z)p £ (t)). 

In the case K = — 1, the coupling to the Poisson equation comprises an attrac- 
tive nonlinearity for the Schrodinger-type equations (|1.2p and hence, global well- 
posedness for general initial data does not hold, cf. [5J[T3]. In this case, the system 
(jl.2p is a generalization (for mixed states) of the semi-relativistic Hartree model 
derived in [2] as the mean field limit for large systems of Bosons with gravitational 
self-interaction. This model has been extensively studied in recent years as it is con- 
sidered to describe the dynamics of so-called Boson stars [H El [T3] • In the case 
k = 1, the Poisson interaction is repulsive and typically models electron-electron 
self-interactions (in the Hartree-approximation). The system (ll.2j) therefore allows 
to study relativistic corrections to the usual Hartree model of many-body electron 
system as it is needed for example in the case of heavy atoms, cf. [14] . 

In the present work we are interested in rigorously establishing the classical limit 
of fll.2|) as e — > 0. To this end we shall heavily rely on the Wigner transformed 
picture of quantum mechanics [21) . To this end, we define the (e-scaled) Wigner 
transform of a given density matrix kernel p e {t, x, y) as in [16[ I15j 

1 f E ( M _ , ey ey N 
(2tt) 3 

From this definition we easily infer 

C 1 - 3 ) ll/ £ (*> ') Ollz^ItgxRf) = (27T£)~ 3 ||p e (t, •, -)IIl2(R3xR3)- 

Since f e (t,x,£,) is real- valued it can be seen as a quantum mechanical analog of 
a classical phase space distribution. However, the Wigner function in general also 
takes negative values and therefore does not allow for a probabilistic interpretation. 
Nevertheless, Wigner functions have proved to be a highly successful tool in the 
rigorous mathematical derivation of classical limits, see, e.g. [H l8l [91, IT5 [ \TE \ H5 1 [20] 
for various analytical results and [17] for a numerical study. In particular, under 
appropriate uniform bounds on p e (see Theorem 1 below), it is known that f e has 
accumulation points / = f (t, x,£) as e — > (in a suitable topology), which are 



f[f{t)] = f(t, x, := _ / 3 f (t,x+^,x- f ) e-^y dy, 
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positive Borel measures on the phase space. The limiting measures can then be 
identified as a distributional solution of the corresponding limiting evolutionary 
system. In our case we expect it to be a solution of the relativistic Vlasov-Poisson 
system, i.e. 

I d t f + ■ j • V x / - V X V ■ Vtf = 0, 
(1.4) VWTT U 

[ - kAV = n(t,x), 
where the limiting particle density n(t, x) is obtained as 

n{t,x) = / f(t,x,£)d£. 



The system (|1.4p is in itself an intensively studied model, in particular in the 
gravitational case k = — 1, see e.g. [TUl [TT1 [T2] and the references given therein. The 
present work thus provides a rigorous connection between the quantum mechanical 
model (|l.ip . or equivalentlv (|1.2p . and its classical analog (|1.4|) . The novelty of 
the work lies in the fact that we have to establish suitable estimates based on the 
relativistic kinetic energy (instead of the usual one — %-A) in order to be able to 
pass to the limit in the nonlinear potential V e (t,x). The main result of this paper 
is as follows. 

Theorem 1. Let Eq < 1 and assume that the initial data p$ E (5i(L 2 (R 3 )) is a 
density matrix operator, such that 

(A) sup (tr^ + ^tr(^) 2 + trv / -e 2 A + l^J < oo. 



ee(0,e o ] V e3 

In the gravitational case k = — 1 we additionally assume 
(B) ^tr(p-f< ( 



trpg 

where C* > is a fixed e -independent constant to be computed in the following. 
Then there exists a unique mild solution p e G C([0, oo); ©i(L 2 (R 3 )) of the relativis- 
tic Hartree system (II. ip and its Wigner transform f £ [p e ] E C([0, oo); L 2 (P 3 
converges, up to extraction of sub-sequences, 

f[p e ] £ -^f inC([0,T};S'(Rl xl|)-w*), 

for anyr < oo, where f E C([0, r]; M + {R 3 X x M 3 )) n L°°([0, r]; L 1 n L 2 (K 3 x M 3 )) is 
a distributional solution of the relativistic Vlasov-Poisson system (|1.4p . Here M + 
is the space of positive bounded Borel measures on phase-space, equipped with the 
weak-* topology. 



The above given theorem shows that distributional solutions to (|1.4|) can indeed 
be interpreted as the classical limit (on any compact time-interval) of solutions to 
(II. 2p . or equivalently (|l.ip . To get more insight on Assumption we again use 
the decomposition 

to obtain 

trpg=X)A|, tr(^) 2 = £(A^) 2 . 

We remark that trpg is the total charge or mass of the particle system. Thus 
Assumption (|A|) implies a uniform (in e) bound on the total charge/mass and 
simultaneously requires ^j^i^j) 2 — £ 3 , as e — * 0. (For a construction of a density 
matrix pg, which satisfies this requirements for tpj E i/ 1 (R 3 ) we refer to 16J.) 
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This condition is to be expected from earlier papers [T^l [H] and prevents us from 
establishing our result in the case of pure initial states, i.e. p§ = ^j(x)ipj(y), or 
even finite combinations of pure states. Roughly speaking, the requirement of a 
totally mixed initial state is needed to ensure that the limiting particle density 
n(t, x) is not too singular and thus can be successfully convolved with the Poisson 
kernel oc l/|x|. This is not clear a-priori as the limiting f(t, x,£), and thus also the 
limiting density n(t,x), in general is only a measure. In the case of a fully mixed 
state though we gain a bit more regularity, which is needed in passing to the limit 
e -> 0. Note that for any Hilbert-Schmidt p e <E 6 2 (X 2 (K 3 )), it holds 

IIP £ le 2 (L 2 (R 3 )) = t r {p e ) 2 — \\p £ (', ")IIl2( R 3 xR 3) < oo , 

cf. [T5]. Having in mind the scaling property Q1.3|) . Assumption ((TfJ therefore 
implies for the initial Wigner function ||/ e (0, •, -)IIl 2 < °°i uniformly in e. This 
property is shown to be conserved by the time-evolution below and thus, yields 
an important uniform bound on X, £). Assumption fBj). needed in the grav- 
itational case, then additionally requires that the L 2 -norm of the initial Wigner 
function or the total charge/mass is sufficiently small (and not only bounded), cf. 
Remark 12.51 below. 

The rest of the paper is now organized as follows: In Section ([2]) we collect some 
preliminary results needed for the proof of our main theorem which is given in 
Section [3] We then close this section with some final remarks on possible general- 
izations. 



2. Preliminary results 

Let us start with a technical lemma, that nevertheless turns out to be crucial for 
the proof of the main Theorem. 

Lemma 2.1. Let p(x, y) be the kernel of a positive self-adjoint trace class operator 
p G 6i(L 2 (R 3 )) and n(x) = p(x,x) the corresponding density. Then for p € [l,oo) 
the following estimate holds 

(2-1) IMIi,(R») < C P IHIe^L^)) (tr Mp) 1 - 9 , 

where 

= 4p- 3 = p 

q ~3p-2 1 ~ 4p-3' 

Proof. For the proof we proceed analogously to [15J. Since the case p = 1 is 
immediate we consider only p > 1. We start by considering the operator | V| — (3n a 
where j3 > and a > are some constants to be chosen later on. Let p\ < p2 < 
. . . be the negative eigenvalues of the operator |V| — [3n a and let {ifk}ken be a 
finite or countable collection of the corresponding ortho-normed eigenvectors. We 
consequently obtain 

tr(|V|-/3n Q )p>^((|V|-/3n Q )p^,^>= E M^fcA,, 

keN j.keti 

where pjk = L 3 ipj^Pk- Since pk = — | , V € N, we can rearrange this inequality 
in the following form 

tr f3n a p = f3 f n a+1 dx< tr|V|p+ ^ \ Pjk \ 2 \pk\\ r 
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Now, using Holder's inequality, and the fact that Y^j \pjk\ 2 = 1, as well as J2 k \pjk\ 2 < 
1, we obtain 



(3 / n a+1 dx < tr|V|p+ ||p|| 



S P (L 2 (R 3 )) 



< te\V\p+\\p\\ 6p{L 2 {R3)) C p 

Here the second inequality is obtained from Theorem 2.1 in [3], which states that 
for all 5 > 




~z1 Js. d 



where p,j are the negative eigenvalues of the operator | V| — W. Now with the choice 
a = (p — l)/(3p — 2), setting a + 1 = q and 



Q_( iSv 



1/p \P'/3 



\2C llplle p (L 2 (R 3 )) 
it is easy to conclude 

i/g 



□ 



Remark 2.2. The proof can easily be generalized to the d-dimensional case, where 
one finds 

= d{p -l)+p = P 

q d(p-l) + l' d(p-l)+p' 

Next, we state a local- in-time existence result for the Hartree system (|1.2[) . 



Lemma 2.3. Lei pg 6e smc/i i/iai Assumption (|A"|) holds. Then there exists a T > 
and unique mild solution p e e C([0, T); Si(L 2 (R 3 )) to ([O]) . or equivalently (fTTTj) , 
which satisfies the following conservation laws: 

(2.2) trp e (i)=tr Po \ tr(p £ (t)) 2 = tr(^) 2 . 

Proof. The existence of a unique solution can be obtained in a straightforward way 
from the results given in [5] (where the case of density matrices with finite rank is 
treated). Also, the conservation laws (|2.2[) are established there. □ 

Note that the conservation laws ()2.2[) together with Assumption j7[J directly 
imply that, for all e E (0, eq], it holds 

IKWIIl1(R3) = 11^(0)1^1^3) < 00, 

and 

II/ E (*)IIl 2 (R3xR 3 ) = II/ £ (°)IIl 2 (R3 xR|) < °°> 

for all t £ [0,T). 

Lemma 2.4. The solution stated above exists globally in-time, i.e. T = oo, and 
we additionally have conservation of energy: £%At) — £|j(0), provided that one of 
the following conditions is satisfied: 

i) K = 1, 

ii) K = — 1 and Assumption fB]) is satisfied. 
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In particular we have the following uniform bound 
(2.3) £un(*)<C, Vte[0,oo). 

Proof. Formally, the conservation of energy can be obtained as in [5] . In order to 
obtain the uniform a-priori bound (|2.3|) we shall use (|2.1|) . with p = 2, which yields 



||n e (t)|| LB/ .<C 2 ||^(t,.,OHi / s 8 (tr|V|p e (t)) 8 / 6 

< C 2 e- 3 ^\\f(0, -)\\% B (tr y/-e*A + l p s {t) 



due to the second conservation law in (|2.2p and the fact that e|V| < e 2 A + 1, 
as can be seen on the level of the Fourier-symbols. On the other hand, recalling 
the Sobolev inequality 

n'{t*W,v) Q >||n . (f)|| a 

;R3 \X-y\ 

and interpolating the density via ||^ £ (i)||^6/5 < ||^ 6 (i)||^i 6 || nE (*)|| r,B/4) we obtain 



// 

J Jr 3 



" (/ '"" (l v) dxdy<C£^ n (t), VtG[0,oo), 



y\ 



where C = C s C^ 3 (tr/3g) 1 ^ 3 ||/ e (0, •, • ) 1 1 i s £-independent by Assumption l[A"|). On 
the one hand, this shows that the initial energy is indeed well defined and uniformly 
bounded in e for both k = ±1. Moreover, we immediately conclude the uniform 
boundedness of the kinetic energy in the case k = 1 and hence, global-in-time 
existence of solutions. In the case k = — 1 we only get, from energy conservation, 
that 

Uniform boundedness of £^ in (t), and hence, global- in-time existence, can therefore 
be concluded if C < 8ir, which holds true under Assumption [Bl □ 

Remark 2.5. In comparison to the proof given in [13] we needed to invoke slightly 
different arguments in order to obtain the uniform (w.r.t. e) bound (12. 3p . The 
requirement C < 8tt directly leads to Assumption jB]) with C* = ^j, s . Note 
however, that in our mixed-state case we can no longer characterize the condition 
C < 8ir by the solution of the corresponding single-state ground state problem as 
it is done in [T5] . 

3. Proof of the Main Result 
In this paper we shall use the following definition for the Fourier transform 

{T^){n) ee 0( V ) = j 3 viOe-^dt. 

Moreover we denote by S the Schwartz space of rapidly decaying functions. 
Proof. The proof of the Theorem [T] consists of three steps: 

Step 1. We first note that due to Assumption and the conservation laws 
(12. 2|) . the Wigner function f e (t) is uniformly bounded in L?(M? X x K|) and in 
S'(M 8 x R|), where the later follows from the results in [3]. Thus, for every fixed 
t G [0, oo), it holds (up to extraction of sub-sequences): f e {t) — ^ fit), as e — > 0, in 
- w* and in L 2 (M.^. x R?) — w. Moreover it has been shown in [51 [T5] 
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that the limit / £ ./V1 + (R 3 x Rjp, i.e. a positive phase-space measure. 



Step 2. Next we shall prove the time-equicontinuity of f £ (t,x,£,). To this end, 
we have to show that d t f £ is bounded in £°°((0, t), S'(J$? x x R|)), for any r < 
oo. We consequently consider the (phase space) weak formulation of the Wigner 
transformed evolution equation, i.e. 



(3.1) 



{d t f £ 1 <j ) ) = {T £ f £ ,<p) + {Q £ [v £ ]f £ 



where <f> £ S(R X x R|) and we write ( • , • ) for the corresponding duality bracket 
on S. Above, T £ and O e [U e ] are pseudo-differential operators, corresponding to the 
Wigner transformed kinetic and potential energy terms. Explicitly, T £ f £ is given 
by 

i 



r s f(t,x,o = 



(2tt) s 



where 



7 e (z/,£) 



2£ ^ ™a 



On the other hand, the nonlinear potential V £ (t, x) enters via the well known 
operator, cf. [IBlfTC] 



e £ [v £ }f(t, x ,o 



(2tt) 3 



whose symbol is given by 
5 £ (t,x,r,) :=±(y e (t,x + el)-V e (t,x-elj) = rj ■ f' 



1/2 



/2 



V X V £ (t,x + ssr}) ds. 



In order to estimate the first term on the r.h.s. of (|3.ip . we use an inverse Fourier 
transform w.r.t. x and Planchcrcl's theorem to obtain 



l(r £ / £ ,0)| 



< 



( J ^-"(J^VvMy, 0-7 e (Z,y) dy)f E (t, v ,Odr) d£ 
f E (t)\\L*( I f |(^yV^)(i/,0| a l7 £ (e,y)| 2 dyde) 



1/2 



< ll/'WIWIv^ 



lL 2 ' 



where the second inequality follows from the fact that |7 e | < 1, for all e £ (0, £o]- 
Since f £ (t, x, £) is uniformly bounded in L 2 (R 3 X x R|) for all t £ [0, oo), due to 
and Assumption |(A[). we obtain the desired uniform bound for |(r e / e ,0)|. Next, 
we consider the nonlinear term on the r.h.s. of (|3.ip . After a Fourier transform 
w.r.t. £, we need to estimate 

rl/2 



|(e £ [n/ £ ,0)l 



-1/2 



-3/2, 



Using the generalized Young inequality we have 

iiv,u £ (t)ii L2(R3) < cljv^r 1 

and by interpolation between L 1 (R 3 ) and L 5 / 4 (R 3 ), one obtains 

\\n £ (t)\\ L6/5 < lln^tJII^^In^tJII^^M^W)^!!/^.,^!^ 3 ^^. 



Due to the conservation laws (|2.2p and the uniform bound on the kinetic energy 
(|2.3p we therefore find that || ^ x V £ {t)\\ L 2 is uniformly bounded as e — > 0, for all 
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t £ [0,oo). Thus, we can estimate 

"1/2 2 ,1/2 



W[V £ ]r,<P)\<\\f{t)\\J II l#^)| 2 M 2 / V x V%t,x + esrj)ds 2 dxd V )' 

v JJl6 J -1/2 

< \\f £ (t)\\ L2 \\V x V s (t)\\ L J f sup \4>(y - eswtflnfdri) 1 

v Jm yew 3 ' 



,V2 

>(2/-£rs77,7?)|"|77|'dry, 

/R3 yGR 3 

< ||/ e (*)|| ia l|V,y £ (t)|| i2 || sup \4>(y, V )\\ V \ IU2 < ex., 



uniformly in e. Thus, d t f £ is bounded in L°°((Q,t),S'(R x x R|)), uniformly in e, 
and we conclude that (again up to extraction of sub-sequences) 

f^f in C((0,t);S'(R 3 x x Rf) - «*), 

and in i°°((0, t); i 2 (M|J x R|)) - w * . 

Together with the results of Step -? this implies 

/ e C((0, t);M+(M.I x R|) - to*) n £°°((0, t);L 1 n £ 2 (R^ x R|)). 

Step 5. It remains to identify the limiting f(t,x,!;), by passing to the limit in 
the (full) weak formulation of the Wigner transformed evolution equation, i.e. 

(3.2)^ 

( - d t a(t)cj)(x, - a(t)r e 0(x, 0) f £ (t, x, dx d£ dt 



+ 



( 7(i)6 £ [n^,0/ £ (i,a ; ,0da ; d^di= /Y a(0)<Kz, 0/o M d*d£, 



where <r € Cjf (R^) and G S(RJ x R|). First it is easily seen that 
Urn sup |(/M^)-(/,V^- Jj— )|=0, 

since / e / in L°°((0, r]; £ 2 (Rj x R|)) - w*, as shown in Step 2, and the conver- 
gence of the rest of the duality bracket is strong in C([0, T]; i 2 (M^ x Rjp). Con- 
cerning the nonlinear term, we have to show 

(3.3) lim sup |(f,e E r]^)-(/,V^.V^)|=0. 

e ^°te[o,T] 

Similarly to what is done above, after invoking a Fourier transform w.r.t. £, it is 
sufficient to prove that 

lim ($(x, rj) I 77 • V x V £ (t, x + esr))ds) = <j>(x, 77)77 • V x V(t, x), 

in C([0, T|; L 2 (R^ x R|)) strongly. In order to do so we denote the compact support 
of <j> by :C Q, x x £l v and write 

/ 1/2 ( IL ^If^ccV'it, x + ssr]) - V x V(t, x)\ 2 dx&n) ^ ds 
< sup \4>(x,r])\r]\\ I (II \V x V £ (t, x + esrj) - V x V(t, x + es?7)| 2 da; d?7 

x.-qeil J -1/2 V J JO. 

+ J J \W x V(t,x + esT]) - V :E F(t,x)| 2 da;d?7 ) ^ds 
= : sup \4>{x, V )\r,\\ (I £ + J 6 ) 1 ^ 2 ds. 

i,i|G!i J -1/2 
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First, we consider the term J £ : Using Plancherel's Theorem we can write 
J £ = // \V x V(t,x + earj)- V x V(t,x)\ 2 dxdr) 

\(^ z V x V)(t,z)\ 2 \e iss ^ -l\ 2 dxd V , 

from which we conclude that J £ — > 0, as e — > 0, by dominant convergence. In order 
to treat the term I £ we take into account that V s solves the Poisson equation 

-kAV £ = n £ (t,x), k = ±1 

with n £ (t) G L 5 / 4 (R 3 ) D L 1 (IR 3 ), uniformly in e. Due to the regularizing property 
of the Poisson equation we have V X V £ G H^' 5 (R 3 ) uniformly in e, which is 
compactly embedded in L[ 2 oc (lR 3 ). We therefore infer that, for all t G [0, r], it holds 

(3.4) V x V £ (t)^V x V(t) in L 2 (0), 

where O — {z := x + esrj G M 3 | (x,r)) G (SI 3 x SI 3 )}. We finally recall the result 
of Theorem III. 2 in p~5] which (up to extraction of subsequences) states that, as 
e -► 0, 

n e (t,x) ^ n(t,x) = f(t,x,£)d£, 

JR 3 

in C([0, T]; A4 + (K 3 ) - to*). We therefore conclude J E — ► 0, as e — > 0, which con- 
sequently implies (|3.3[) . In summary we have shown that the Wigner transformed 
evolution equation (|3.2[) converges weakly to 

( - d t a{t)cj>{x, - -jJ— ■ V x 4(x, £)a(tj) f(t, x, f) do: d£ dt 
v V / I?l+ 1 7 

+ / // <r(t)V x V-Vt<t>(x,Z)f(t,x,Odxd£dt= [[ a(0)ci)(x,Ofo(x,t;)dxdZ, 

JO J JRO J JR6 

i.e. the relativistic Vlasov-Poisson system (in distributional sense). This finishes 
the proof. □ 

Remark 3.1. It is straightforward to generalize our result for Yukawa-type inter- 
actions, i.e. 

g— A|£C 

V £ (t,x) = n — -*n £ (t,x), A > 0, 

Air\x\ 

which arc used in, e.g. 13J . Indeed it is possible to prove our theorem for an even 
more general class of self-consistent interaction potentials V £ = G * n £ , if the kernel 
G = G(x) satisfies appropriate regularity conditions, cf. [TJ Q~5] for a discussion 
in the case of non-relativistic Hartree systems. Finally, we remark that one can 
also treat the case of semi-relativistic Hartree-Fock systems, cf. [5 , by combining 
straightforwardly the results present here with those given in [7] for non-relativistic 
Hartree-Fock systems. To this end, one should note that the so-called exchange- 
term vanishes as e — > 0. Hence, one again recovers the relativistic Vlasov-Poisson 
system (|1.4[) in the limit. 
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